Extracting signatures of quantum criticality 
in the finite-temperature behavior of many-body systems 
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We face the problem of detecting and featuring footprints of quantum criticality in the finite- 
temperature behavior of quantum many-body systems. Our strategy is that of comparing the phase 
diagram of a system displaying a T = quantum phase transition with that of its classical limit, in 
order to single out the genuinely quantum effects. To this aim, we consider the one-dimensional Ising 
model in a transverse field: while the quantum 5=1/2 Ising chain is exactly solvable and extensively 
studied, results for the classical limit (S — > oo) of such model are lacking, and we supply them 
here. They are obtained numerically, via the Transfer-matrix method, and their asymptotic low- 
temperature behavior is also derived analytically by self-consistent spin-wave theory. We draw the 
classical phase-diagram according to the same procedure followed in the quantum analysis, and the 
two phase diagrams are found unexpectedly similar: Three regimes are detected also in the classical 
case, each characterized by a functional dependence of the correlation length on temperature and 
field analogous to that of the quantum model. What discriminates the classical from the quantum 
case are the different values of the exponents entering such dependencies, a consequence of the 
different nature of zero-temperature quantum fluctuations with respect to thermal ones. 



I. INTRODUCTION 

One of the most fascinating aspects of many-body sys- 
tems is the possible occurrence of a phase transition, ei- 
ther at finite or at zero temperature. This latter case is 
generally referred to as a genuine quantum phase tran- 
sition (QPT), by this meaning that it is exclusively ob- 
served in quantum many-body systems. Zero-point quan- 
tum fluctuations are recognized as the fundamental in- 
gredient of a QPT, in the same sense as thermal fluc- 
tuations are in ordinary finite-temperature phase transi- 
tions: Whenever fluctuations are frozen, as in the T = 
classical case, no phase transition may possibly occur. 

In the literature, a QPT is also commonly said to occur 
when, for a given value of one of the Hamiltonian param- 
eters, the ground state of the model qualitatively changes 
its structure. This definition is not quite rigorous, firstly 
because it labels QPT any change in the universality class 
of the model, as well as mean-field phenomena such as 
saturation, secondly because it may paradoxically be ex- 
tended to classical systems: As a matter of fact, qual- 
itative changes in the structure of the minimum-energy 
configuration may well be observed, for a given value of 
some Hamiltonian parameter, also in classical systems at 
T = 0, despite fluctuations being frozen. In this frame- 
work, therefore, the comparison between the behavior of 
a quantum system displaying a QPT and that of its clas- 
sical limit becomes meaningful even at zero temperature, 
as shown in the next Section. 

When temperature is switched on, the very same 
definition of QPT looses its meaning; however, it is 
well established that a QPT induces a peculiar finite- 



temperature behavior which is qualitatively described by 
the best known phase diagram introduced in Ref. [l], and 
reported in Fig. [3 where three different regimes appear 
(renormalized classical, quantum critical, and quantum 
disordered), separated by crossover lines thoroughly dis- 
cussed in Ref. The relevance of this phase diagram 
is mostly due to its suggesting that signatures of a gen- 
uine quantum critical behavior may survive also at fi- 
nite temperature, a fact that opens the possibility to 
observe them experimentally. Moreover, a renewed in- 
terest has arisen since entanglement properties have en- 
tered the physics of many-body systems, and questions 
like "how resistant to thermal noise are certain quantum 
properties?" became essential in order to test possible 
realizations of quantum devices. 

Despite the above phase diagram being considered as 
strictly peculiar to quantum systems, especially as far 
as the quantum critical and quantum disordered regimes 
are concerned, its structure results from the interplay be- 
tween thermal (classical) and quantum fluctuations: A 
precise analysis of the role played by these two compo- 
nents is therefore necessary in order to ascertain whether 
the latter play an essential role or not, and to distill gen- 
uine footprints of quantum criticality to be experimen- 
tally looked for. To this end, knowing the behavior of 
the quantum system at finite temperature is not enough, 
and a careful comparison with the corresponding classical 
limit is necessary. The quite unexpected lack of results 
for the classical limit of quantum models displaying a 
QPT has made such comparison unavailable until now. 

This paper is aimed at filling this gap: We consider one 
of the paradigmatic models displaying a QPT, namely 
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the one-dimensional quantum Ising model in a trans- 
verse field (QIF), and compare its behavior with that 
of its classical limit, namely the classical Ising model in 
a transverse field (CIF), which we study numerically, via 
the Transfer-matrix method, and analytically, via self- 
consistent spin-wave theory. The result is unexpected: A 
finite temperature phase diagram is disclosed also for the 
CIF, and it has the same structure of that for the QIF. 
In full analogy with the quantum case, we identify three 
different regimes on the basis of the field and tempera- 
ture dependence of the correlation length; moreover, we 
show that the algebraic behavior, which was thought to 
specifically characterize the quantum critical and quan- 
tum disordered regimes, does in fact show up also in the 
classical system, though with different exponents. This 
result jeopardizes the experimental renderings based on 
the statement that the observation of an algebraic depen- 
dence of the correlation length implies the occurrence of 
a quantum critical or disordered regime 3 ' 4 ' 5 ' 6 . In fact, 
an accurate analysis of the exponents is here shown to 
be necessary in order to discriminate genuine quantum 
effects. 

The structure of the paper is as follows: in Sec. ITT1 we 
introduce the model and discuss its zero-temperature be- 
havior. In Sec. IIII) we first summarize the known results 
for the QIF, and then present our results for its classi- 
cal limit, the CIF: numerical data for the magnetization 
and the susceptibility in the field direction, and for the 
specific heat are shown and discussed. The analysis of 
the field and temperature dependence of the correlation 
length is considered in Sec. IIV1 where the classical phase 
diagram is finally obtained and compared with that for 
the QIF. Conclusions are drawn in Sec. El 



II. THE ISING MODEL IN A TRANSVERSE 
FIELD: T = 

A. The quantum model 
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FIG. 1: Zero temperature limit of m x , m 2 , and £ x of the 
classical (solid lines) and of the quantum (dashed lines) Ising 
model in a transverse field. 



and 



(H c - Hf for H -> H c 



Cx ~ (H - H c )- V for H^H+ 



(2) 



(3) 



with the exponents j3 = 1/8 and v = 1, as in the 
finite-temperature phase transition of the classical two- 
dimensional Ising model. 10 Together with the other usual 
exponents they obey typical scaling relations which, due 
to the intrinsically dynamical nature of quantum fluctu- 
ations, entail^ the dynamical critical exponent z. On the 
other hand, the magnetization along the field direction, 
m z = (Sf)/S, is an analytic function of the field and 
changes its curvature at H c , where the uniform suscep- 
tibility consequently displays a maximum. These T = 
behaviors are reported as dashed lines in Fig. [TJ 



One of the best known examples of a many-body sys- 
tem displaying a QPT^l is the one-dimensional quantum 
Ising model in a transverse field (QIF), whose Hamilto- 
nian reads 



n 
j 



i+l 



HSf 



(1) 



where i runs over the sites of an infinite chain, and Si are 
£=1/2 spin operators, J is the exchange energy constant 
and H is the transverse field in units of J; this model is 
exactly solvable^ by means of a Jordan- Wigner trans- 
formation to Fermi operators, and displays a QPT at 
T = and H = H c = 1/2: The discrete Sf -> -Sf global 
symmetry of the Hamiltonian is spontaneously broken 
for H < H c , where the order parameter m x — (Sf)/S 
becomes nonzero; long-range order sets in at the crit- 
ical point, as testified by the divergence of the order- 
parameter correlation length. In particular, it is 



B. The classical model 

The classical limit of any given quantum system is 
unique, although the converse is not true. For a spin 
system the limit h— >0 must be taken while keeping finite 
the spin angular momenta hSi. 

From Eq. (fTJ) one thus obtains the Hamiltonian of the 
classical Ising model in a transverse field (CIF), 



H 

Jc 



N 



V l+ ir2ft<). (4) 

i=l 

where Si — (sf,s^,sf) are classical spins, i.e., three- 
dimensional vectors of fixed length |s;| 2 = 1, whileii 



J c = lim JS , 

S — >oc 



2h = lim — 

R->0 S 
S — >oc 



3 



are the classical exchange interaction and reduced field, 
respectively. Taking the exchange interaction as the en- 
ergy unit, the dimensionless temperature will be denoted 
with t = T/J c , so that all thermodynamic quantities de- 
pend upon the pair (ft, t). Periodic boundary conditions 
are assumed (sn+i — s\) and the thermodynamic limit 
(N — > oo) will be considered. 

It is important to distinguish between the CIF and 
what is often called the classical Ising model, which is 
obtained from Eq. {T]) with H = 0, by regarding the quan- 
tum operators Sf as classical discrete variables taking the 
values ±1/2. Results for the CIF are lacking, apart from 
the zero-field cas o 12 i 13 ' 14 , probably because, at variance 
with the QIF, the CIF does not allow for an exact solu- 
tion. 

Writing the classical spin variables in terms of polar 
angles as 

Si = ( sin 6*j, cos 0j sint^, cos 0j cos ipA , (5) 
the Hamiltonian (|U) is expressed as 

— = — ( sin 9i sin 9 i+ i + 2ft cos 0i cos i/jj) ; (6) 

i 

its minimum corresponds to a translation-invariant con- 
figuration { Si = ( sin 9 m , 0, cos 6 m ) } , with 



± cos 1 h for h < 1 







for ft > 1 



The minimum energy per spin is 
u(h,0) -- 



-1-h 2 for h < 1 
-2ft for ft > 1 



(7) 



(8) 



and shows a singularity at ft = 1 . The magnetization m z 
is proportional to the field for ft < 1 and saturates for 
ft>l, 



m 2 (ft,0) 



ft for ft < 1 
1 for ft > 1 



(9) 



so that the corresponding susceptibility is discontinuous 
at ft = 1 



Xz(h,0) = d h m z 



1 for ft < 1 
for ft > 1 



(10) 



Finally, the behavior of the magnetization along the ex- 
change, 



m x (ft,0) = 



±Vl - ft 2 for ft < 1 







for ft > 1 



(11) 



reflects the fact that for ft < 1 the minimum is twofold. 

It may sound odd but, as seen in Fig. 1, the CIF dis- 
plays a zero-temperature behavior which is analogous to 
that observed in the QIF, even if no fluctuations are 




FIG. 2: Second derivative of the specific heath of the QIC, 
d 2 H c(H,T). 



present at t = 0. Eq. (fTTj) shows indeed that a criti- 
cal field, h c — 1, separates a symmetry-broken minimum- 
energy configuration with m x ^ from one with m x = 0. 

An even closer analogy is found if one considers that 
from the low-temperature expression derived in Ap- 
pendix [B] one can obtain the exact zero-i limit of the 
correlation length, £2: (ft, 0) = 1/y/h—l. The counterparts 
of both Eqs. (J2|) and ((3]) are then available, and read 

hf for ft ->■ ft" , 



and 



(ftc 



Zx~(h-h c )- v for ft ^ 



(12) 



(13) 



with Gaussian critical exponents (3 
compared with those for the QIF, (3 ■ 



-- v = 1/2, to be 
1/8 and v = 1. 



III. THE ISING MODEL IN A TRANSVERSE 
FIELD : T > 

A. T > : The quantum model 

The field- and temperature dependence of the specific 
heat c(H, T) and of the susceptibility Xz {H, T) can be 
easily obtained from the analytic results of Ref. [^. The 
most prominent feature is the occurrence of maxima of 
both quantities in the H—T plane. Indeed, the quantum 
specific heat just shows the behavior of a free Fermi gas 
with dispersion ojk ~ A + fc 2 in the neighborhood of a 
vanishing gap A ~ \H — H c \, which appears, for fixed 
low-T, as two symmetric peaks at linearly displaced po- 
sitions \H — H c \ cx T. This feature is made evident in 
Fig. [21 by the density plot of the second derivative of 
c(H, T) with respect to H . The positions of the max- 
ima draw two symmetric lines in the H — T plane which 
coincide with those obtained^ from the analysis of en- 
tanglement properties. 



B. T > : The classical model 

The thermodynamic behavior of the CIF for h < h c 
is essentially determined by the energy landscape of 
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the model. In a mean- field approach, i.e., setting 
{9i = 9, (fi = 0}, this is described by the double- well en- 
ergy profile e(0) = — sin 2 9 — 2hcos8, with minima in 
9 m = ± cos -1 h and the barrier top at 6 — 0, the barrier 
energy being 5e = (l — h) 2 = (h c — h) 2 . The two wells 
correspond to the Ising configurations. 

Domain- wall excitations connecting the two Ising con- 
figurations can appear on the chain, the energy of a do- 
main wall being e w = 2(1 — h 2 ). A simple statistical argu- 
ment gives a finite number of domain walls at any finite 
temperature, n w ~iV/(l + e 6 ™/*), so that when temper- 
ature is switched on the ordered state is destroyed (i.e., 
m x = for t > 0) by these excitations, which rule the low- 
temperature thermodynamics in what we will hereafter 
call the Ising regime. 

When the temperature reaches the order of the bar- 
rier energy, t~ {h c — h) 2 , thermally activated transitions 
between the wells can occur and the Ising regime breaks 
down. The dependence of the barrier height on the field 
is responsible for the fact that as h — *■ h~ this regime 
gets confined into the narrow interval < t < (h c — h) 2 . 
The above mechanism, that we call thermal hopping, is 
at the hearth of the phenomenology of the model below 
h c , and it already suggests the occurrence of a crossover 
from an Ising-like behavior towards a critical one, ruled 
by an effectively flat energy-landscape. 

On this basis, let us discuss the finite-temperature data 
obtained^ by the Transfer-matrix method, briefly de- 
scribed in Appendix |A"1 

The temperature behavior of the magnetization along 
the field, m z (h,t), is shown in Fig. [3l and clearly is re- 
lated to that of the mean-square fluctuations of the Ising 
order parameter, ((sf ) 2 ); indeed, the initial increase with 
temperature corresponds to the reduction of ((sf ) 2 ) due 
to the fact that the probability distribution of sf , initially 
frozen in the bottom of one of the two wells (±vl — h 2 ), 
extends more likely towards the barrier at sf = 0. The 
further decay of m z is due to the isotropic spin fluctu- 
ations occurring after thermal hopping has taken place. 
On the other hand, for h>h c , m z simply decreases from 
its t = saturation value. 

The susceptibility x.z{h, t) is shown in Fig.|4j the zero- 
field result agrees with that derived in Ref. [lj and shows 
a broad maximum in temperature, at t ~ 0.37. Upon ris- 
ing the field, such maximum is squeezed towards lower 
temperatures, meanwhile getting sharper. At the criti- 
cal field, the maximum disappears and the susceptibil- 
ity is a monotonic function of temperature for whatever 
h>h c . The zero-t limiting value of Xz(h,t) is given by 
Eq. (|10p. The overall behavior of Xz{h, t) in the h-t plane 
is evidently characterized by the occurrence of the above 
described maxima for h <h c . 

The specific heat c(h, t), shown in Fig. is also char- 
acterized by the occurrence of maxima for h < h c , which 
disappear above the critical field; noticeably they fall into 
almost the same positions as those observed in Xx{h,t), 
as it appears in Fig. [51 

The maxima observed both in the susceptibility and in 
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FIG. 3: Magnetization along the field direction, m z (h,t) vs 
temperature, for selected field values. The 'critical' value is 
h c = l. 
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FIG. 4: Magnetic susceptibility X^{h,t). 




FIG. 5: Specific heat c(h,t). 
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FIG. 6: Region h < h c of the phase diagram of the CIF. Cir- 
cles and squares indicate the position of the maxima of c(h, t) 
and x(h,t), respectively. The dashed curves are obtained fit- 
ting the data with t oc (1 — h) 3 ^ 2 . Inset: log-log plot of the 
same data and curves. 



the specific heat correspond to the onset of thermal hop- 
ping, and their positions indicate the region where ther- 
mal fluctuations overcome the Ising domain-wall energy, 
i.e., the crossover region from the Ising to the critical 
regime. In order to better characterize the correspond- 
ing crossover line we have fitted the maxima positions for 
low t in the h—t plane with the function t oc (h c —h) K , find- 
ing the exponent quite close to the value k = 3/2 derived 
by analytical arguments in the next section. 



IV. THE PHASE DIAGRAM 

A. Quantum model 

The finite-temperature phase diagram for the QIF, 
shown in Fig. [7J essentially features the occurrence of 
three regions, characterized by qualitatively different be- 
havior of physical observables with respect to h and t. 
These regions have been identified^ with the so called 
renormalized- classical (A), quantum- critical (B), and 
quantum- disordered (C) regimes, which have been sin- 
gled out according to the qualitatively different behavior 
of the correlation length £ x , which is 



e u(H c -H)/T 
J 1 — l/l/ 



for H c - H > T (A) 
for \H — H c \ <C T (B) 
for H - i? c > T (C) 



(14) 



where u is a constant and the critical exponents are y — 1 
and v=l. The power law divergence in Eqs. Q3])-(B) 
and -(C) follow from the temperature-dependent scaling 
law for £ x lying at the hearth of the Renormalization 
Group (RG) approach to critical phenomenal. We re- 
mind that the exponent y, ruling the scaling of energy, is 
bound to equal the dynamical critical exponent z by the 
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FIG. 7: Phase diagram of the quantum Ising model in a 
transverse field. 



uncertainty principle. In the same framework, the above 
regimes are shown to be separated by the crossover lines 
Toc\H-H c \. 

This phase-diagram has been extensively discussed 
(see, e.g., Refs. @ and [2) over the last decade: in par- 
ticular, the fact that £ x is independent of the field and 
inversely proportional to T in region B, and independent 
of temperature and inversely proportional to H — H c in 
region C, has always been considered as a signature of 
the genuinely quantum character of the corresponding 
regimes, which have consequently been labeled as "quan- 
tum" (critical and disordered, respectively). 

In fact, as we argue in the remaining part of this article, 
this is not truly the case. 



B. Classical model 



As seen in Sec. IIII B[ hints of the possible occurrence 
of at least two different regimes, separated by a crossover 
region, do already come from the behavior of the suscep- 
tibility and of the specific heat for h < h c . However, in 
order to closely mimic the procedure followed in draw- 
ing the quantum phase diagram, we analyze the field 
and temperature dependence of the correlation length 
£ x (ft, t) of the classical model, heading towards expres- 
sions analogous to those of Eqs. (fT4"]) . The overall behav- 
ior of ^ x {h,t), as from our numerical data, is plotted in 
Figs. E and M 



1. Correlation length for h < h c 

If one considers the ^dependence for h < h c , it clearly 
appears from Fig. [5] that ln^ cx l/t, with a slope that 
increases with the difference h c — h. As for the field de- 
pendence, the logarithm of £ x (h, t) is seen in Fig.[9]to be 
linear in h, with negative slope: this is the more evident 
the smaller the temperature and the angular coefficient 
decreases with increasing t, consistently with Eq. (fT4"|) - 
(A). In fact, from Fig. [10] for h c — h ^> t we see that 
t \n£ x (h,t) w f(t) + a(h c — h), with f(t) weakly depen- 
dent on t, meaning 

£ x (h,t) ~ e ^o-h)/t f or h c -h^>t (a). (15) 
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FIG. 8: Correlation length £ x (h,t) vs 1/t, for selected field 
values. 
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FIG. 9: Correlation length £ x {h,t) vs field, for different tem- 
peratures. The vertical dotted line marks the critical field 
h c = l. 



2. Correlation length for h > h c 



FIG. 10: Classical phase-diagram region a. The function 
t ln^a;(/i, t)] vs h c — h for different fixed temperatures. One 
can see that the curves become straight lines the better the 
condition h c —h^t is satisfied. 




FIG. 11: Classical phase-diagram region b. Longitudinal cor- 
relation length £, x (h, i) vs t for fields close to the critical value 
h c = 1. The log- log plot emphasizes that the slope of the curve 
for h — h c is -1/3. 



From Figs. [8] and [9] we see that for field close and above 
the critical value, the temperature- and field dependence 
of £,x{h, t) becomes much less pronounced with respect to 
those displayed below the critical field, suggesting power- 
law behaviors of the same type observed in the quantum 
regimes B and C. 

In particular, using our numerical data, we can ascer- 
tain that, for \h — h c \ <g; t the correlation length behaves 
as t; x oc i -1 / 3 , as evidenced in Fig. [11] on the other hand, 
for \h — h c \ ^ t, the log- log plot reported in Fig. [T^] em- 
phasizes a power-law field-dependence, £ x oc (h — /i c ) -1 / 2 . 
These behaviors are fully analogous to the quantum ones 
described by Eq. HU)-(B) and -(C), the only difference 
being in the exponents. 

In order to strengthen this result, we develop the 
analytical treatment reported in Appendix [B] which 



makes use of self-consistent spin-wave theory (SSWT) 
for h>h c . The asymptotic behavior of the correlation 
length derived there, valid at low temperature and close 
to criticality, is given by Eq. (|B9[) and yields 

CAh,t)^t-^ 3 for |/i-/i c |«f (b), (16) 

and 

&(M)~ {h-K)- 1 ' 2 for \h-h c \^t (c), (17) 
in full agreement with the analysis of our numerical data. 

3. Crossovers 

Let us now consider the problem of identifying the 
crossover regions between phases where Eqs. (fT5|) . (|16l) . 
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log(h-h) 

FIG. 12: Classical phase-diagram region c. Log-log plot of 
£ x (h,t) vs h—h c , for different low temperatures. The slope of 
the \ow-t curves is seen to be -1/2. 



and (|17p hold, i.e., between the three classical regimes, 
that we have labeled a, b, and c, respectively. 

As for the ab crossover, a reasonable localization can 
be obtained by relating it with thermal hopping. Ac- 
cording to the mean-field analysis presented in Sec. lIIIBl 
the latter occurs when the temperature overcomes the 
energy-barrier height, i.e., for h c — h ~ t x l 2 . However, 
the mean-field approach neglects correlated fluctuations, 
while to get a correct estimate it is necessary to keep 
nearest-neighbor fluctuations at least within a quadratic 
approximation of the Hamiltonian ^ around one of the 
minima. Setting di = m + £•£ and expanding, one finds 
for the quadratic part 



) + h 2 



k 



2(1 - h 2 cos k) e\ , (18) 



from which the mean square fluctuation of Oi around the 
minimum results 



t 



2N ^ 1 

fe 



t 



cosfc 2\/l — h 4 



(19) 



As soon as these fluctuations reach the size of the width 
of the barrier, |0 m | = sin -1 \/l — h 2 , the Ising excitations 
disappear and the crossover between the exponential- and 
the power-law correlation length is expected. This con- 



dition is fulfilled when 
the critical field h c = l, 

t~(hc- h) 



i.e., for h close to 



3/2 



(ab) 



(20) 



This prediction gives the correct exponent and is also 
fully consistent with that obtained in Sec. lIIIBl bv fitting 
the positions of the maxima observed, in the h — t plane, 
for both the susceptibility and the specific heat. 

It is worth mentioning that thermal hopping is the ul- 
timate cause of the finite-temperature phase transition 



to the ordered state occurring for h < h c in more than 
one dimension: the crossover turns indeed into the sharp 
critical line t c (h), with the critical temperature vanishing 
at the quantum critical point, t c (h c ) = 0. In one dimen- 
sion the onset of ordering is forbidden at finite t and only 
the broad crossover between regimes a and b survives. 

Coming to the be crossover, from Eqs. (fTrj)) and (fT7|) 
one deduces it to occur for 



t~{h-h c f' 2 



(be) 



(21) 



4- Phase diagram 

Eqs. (jTSJ) , (fTrj|) . and (fT7|) . can be summarized as: 
e u(h c -h)/t for h c -h->t (a) 



6* 



t -i/y 
(h - h c 



for | h - 
for h 



h c \ < t 
h c > t 



(b) 
(c) 



(22) 



with y = 3 and v = 1/2. Via the different behavior of the 
correlation length, three different regimes can thus be sin- 
gled out also in the classical phase diagram and, accord- 
ing to our discussion, it makes sense to call them Ising 
(a), critical (b), and disordered (c). Different regimes are 
again separated by crossover lines, which in the classical 
model are described by the relation t ~ \h — h c \ z / 2 . 



C. Classical vs quantum phase diagram 

Eqs. (|2"2"]l clearly show that the CIF is characterized 
by a phase diagram on the h — t plane which is fully 
analogous to the celebrated quantum one described by 
Eqs. (|14[) . In particular, the observed power-law diver- 
gence of £ x in regimes (b) and (c) suggests the scaling 
hypothesis to hold also in the classical case: This is con- 
firmed in Appendix|B]where we obtain an explicit expres- 
sion for £ x (h,t), which turns out to be a homogeneous 
function. 

The nexus between the quantum and the classical case 
can be drawn as follows: Let us introduce the unifying 
parameters ~g and t, defined as g = H — H c ,t = T / J in 
the quantum model, and g = h — h c ,t — t in the classical 
model, and write the equation 



(23) 



which, within RG, rules the scaling of observable quan- 
tities and model parameters in the proximity of critical 
points after a length-scale transformation by a factor b. 
From Eq. follow both the quantum Eqs. ©-(B) 
and -(C) and the classical ones, (|22p-(b) and -(c), as 
well as the crossover lines t cx g uv : The critical exponents 
entering the above expressions, despite getting different 
values in the quantum (y — 1, v = 1) and in the classical 
(y = 3, v = 1/2) case, consistently fulfil the hyperscaling 
relation 



2 — a = v{d + y) , 



(24) 



8 



where 2 — a is the exponent for free-energy density /, 
defined^ by f(g, 0) - \g\ 2 -° (a = in both cases). We 
remind that, in the quantum case, the scaling exponent of 
energy y is unavoidably related to the dynamical critical 
exponent z by the uncertainty relation, y — z. As for the 
classical case, we notice that v takes the typical Gaussian 
value, due to fluctuations freezing as t — ► 0. 



APPENDIX A: TRANSFER-MATRIX FOR THE 

CIF 

We outline here the numerical Transfer-matrix tech- 
nique^ by which we have investigated different static 
thermodynamic quantities of the CIF, focusing the at- 
tention on their low-temperature behavior in the neigh- 
borhood of the point h = h c . 

Using polar coordinates one can map the classical spins 
appearing in the Hamiltonian ((4]) as 



CONCLUSIONS 



In this paper we have compared the finite-temperature 
phase diagram of a quantum model displaying a QPT, 
namely the 5=1/2 Ising chain in a transverse field, with 
that of its classical limit. To this end, we have obtained 
numerical and analytical results for the classical model. 
In particular, we have studied the magnetizations, the 
specific heat, the magnetic susceptibility, and the corre- 
lation length along the exchange direction: all quantities 
have been analyzed below, at, and above the saturation 
field, with the temperature raised from zero up to values 
of the order of the exchange interaction. 

The classical phase diagram emerging from our work 
is fully analogous to that of the quantum model. Three 
regimes are identified: the Ising regime a (h c —h^S>t), 
corresponding to the quantum renormalized classical 
regime, where £ x behaves exponentially with (h c —h)/t; 
the critical regime b (\h—h c \ <Ci), corresponding to the 
quantum critical regime, where £ x behaves algebraically 
with t; the disordered regime c (h—h c ^>t), correspond- 
ing to the quantum disordered regime, where £ x be- 
haves algebraically with h—h c . Two crossover lines, 
t~ \h— h c \ 3 / 2 , separate the regions of the phase diagram 
where the above regimes occur. 

The essential message of this work is that in order to 
discriminate quantum critical effects it is not sufficient 
to observe, say, an algebraic behavior of the correlation 
length with respect to temperature, but a precise deter- 
mination of the exponent is rather due. Our analysis does 
also suggest that the role of genuinely quantum fluctu- 
ations at finite temperature is not as relevant as com- 
monly believed, given the fact that most of the features 
regarded as typical of the quantum model are disclosed 
also in its classical limit, even, and most noticeably, at 
very low temperature. 

Given the very weak model-dependence of the overall 
discussion, we believe that the above conclusions hold in 
general, and not only for the Ising chain in a transverse 
field. 
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= (xi, y/l — x? sinifi, a/1 - x 2 cosifi) 



(Al) 



with Xi G [—1,1] and tpi € [0,27r]. The partition function 
can then be expressed as the trace of the AT-th power of 
an integral kernel K(x, y), 

z = y\ f 1 — r dipi c (xixi+i+2h ^^ z ^ cosvz)/t 



i 2 27T 



= f da K N {x,x) =J2 X e > 
J-i i 

where the kernel is real and symmetric, 

cy/t 



(A2) 



K(x,y) = 



l-z* Jo f 



(A3) 



Io(x) is the modified Bessel function, {A^} = {Ao, Ai, ...} 
are the (positive) eigenvalues of K (say, in decreasing 
order), and {^(a;)} the corresponding (real) eigenfunc- 
tions. The diagonalization of K was performed numeri- 
cally after discretizing the integral with a 5-point Simp- 
son's formula 18 on a mesh of up to 1040 intervals, also 
accounting for the definite parity of the eigenfunctions. 

In the thermodynamic limit the free energy per site is 
a function of the largest eigenvalue only, 

f(h,t) = -t lim -rzkiZ N = -t \n\ (h,t) ; (A4) 

the internal energy u — t 2 d t In Ao and the specific heat 
c = dtu, were obtained by numerical differentiation (5- 
point Lagrange formula^). 

The probability distribution for the variable X — x% IS 



w{x) 



^ 2 {x) 



(A5) 



one can reduce to averages with this probability both 
the expressions for the magnetization along the field 
rn z = (s|) = (t/2)dh In Ao and for the corresponding sus- 
ceptibility Xz = dhm z . On the other hand, the joint prob- 
ability for two sites at a distance r, x = Xi and y — x i+ri 



is 



w r {x,y) = Y 



ipoixtygix^oiy^iiy). (A6) 
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Using this result, one can write the correlation function 
of the spin components in the direction of the exchange 
as 



3 i a i+rl 



E 



(A7) 



as the eigenfunction ipQ (x) is even, for large r the leading 
term is that for I = 1 , 



AiV 



3 i i+r7 



(A8) 



so that the corresponding correlation length is given by 



&(M) = [MAo/Ai)]' 



(A9) 



APPENDIX B: SELF-CONSISTENT SPIN- WAVE 
THEORY 

For h > h c the minimum configuration of the classi- 
cal Hamiltonian Eq. (J6j) corresponds to the saturation 
one, i.e., {8i = 0, ipi — 0} . In order to estimate the low- 
temperature behavior of the correlation length, we use 
self-consistent spin- wave theory (SSWT), i.e., we assume 
a Gaussian distribution po = e~^ Wo in terms of a trial 
quadratic Hamiltonian Ha whose coefficients are self- 
consistcntly determined by requiring the identity of the 
po averages of H and Ho, as well as of their first and 
second derivatives. 

In terms of the relevant Gaussian variances D = (9 2 )q, 
D' = (6i9 i+l )o, E = ((f%) (by symmetry (<#i?j) = 0), the 
SSWT amounts to set 



sin W, sm I 



'i+l 



7i+l 



) — cos( 



'i+l 



)] 



— D 



(1+-D) sinh D' - D' cosh D' 



sinh£>' 



-i±I +cosh£>' BiOi+i) 



and 



cos 9i cos ifi = -[cos(9i-<pi) + cos(9i+(pi)] 



ft 



l+F 



where F = (D + E)/2 and we used the SSWT identity 
cosx ~ e-^/ 2 [l + \{{x 2 ) - x 2 )}. 



The SSWT Hamiltonian is diagonal in Fourier space: 

(Bl) 



^ = -Ne (t) + J2 [(A + B^ k )\8 k \ 2 +C\cp k \ 2 



with Hk = 1 — cos fc and 



.4 
B 
C 



he — e 
e-^coshD' 
he- F , 



D-D' 



(B2) 



while eo(t) collects the uniform contributions. It follows 
that the self-consistent expressions for the variances are 



D 



t 



2N ^ A 

k 

t 



1 



t 



D' = —Y-r 

2N A 

k 



+ Bp k 
cos k 



2 V / A{A + 2B) 



= D- 



Bn k ' 2N ^ A 



-y 

9 AT L^i 



li k 



E = 



t 

2C 



B [i k 
(B3) 



The stability condition for Hq is A > 0, which defines 
a threshold field, 



h (h,t) 



e F -°- D ' < 1 



(B4) 



above which the SSWT is meaningful. Note that the 
present approach can describe the finite-temperature be- 
havior of the system also for h < h c = 1 , because the con- 
figuration density can still be approximated by a Gaus- 
sian centered in {6^ = 0, tfi = 0} as long as thermal fluctu- 
ations are large enough to overcome the barrier between 
the two symmetric minima 9i = ± 9 m of H, as explained 
in Sec. ITTlBl 

From the SSWT Hamiltonian (|B1[) the Fourier trans- 
form of the correlation function is immediately found, 



G x (k) = (\8 k \ 2 ) Q 



t 



2{A + B jtzfc) ' 
and can be used to evaluate the correlation length £ 



(B5) 



F 2 



am 

2G(k) 



k=0 



B 
2A 



,F-D 



coshD' 



h — ho(h, t) 



(B6) 



Therefore, when the field is close to the critical value 
h c = 1, the behavior of the correlation length is given by 
£,x ~ [h — hoQift)] -1 ' 2 . In this region the variance D 
is enhanced (A ~ 0) and it can be easily seen seen that 
D' ~ D and F ~ D/2, so that 



A(h,t) = 1 - h (h,t) 
3f 



D + D' - F 



-D 



[(h-l + A)(h + l + A 
~ ct (g + A)- 1 / 2 , 

with g = h— 1 and the constant c = 3/4 v2- 
Rewriting this equation as 



2 



-1/2 



(B7) 



A 

£2/3 



g A 

t 2/3 + £2/3 



-1/2 



(B8) 



it appears that A/t 2 / 3 = F(x), with x = g/t 2 / 3 and the 



asymptotic behaviors F(0) - 



r-2/3 



F{x— >oo) 



-1/2 



From Eqs. (|B6[) and (|B7|) the leading behavior of the 
correlation length is given by 



Z x ~{ g + A)- l / 2 ~t- x f*F{x) , 



(B9) 
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and it follows that near the critical point it is a homoge- which coincides with Eq. (|23|) with the exponents v=l/2 
neous function, and y — 3. For the crossover line (be), identified by im- 

posing x ~ 1, one finds t ~ g 3 / 2 . 

Ub 2 g,bh)^b-^ x {g,t) , (bio) 
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